home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / dswap.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.6 KB  |  74 lines
  1.       SUBROUTINE DSWAP( N, DX, INCX, DY, INCY )
  2. *
  3. *     interchanges two vectors.
  4. *     uses unrolled loops for increments equal one.
  5. *     jack dongarra, linpack, 3/11/78.
  6. *
  7. *     .. Scalar Arguments ..
  8.       INTEGER           INCX, INCY, N
  9. *     ..
  10. *     .. Array Arguments ..
  11.       DOUBLE PRECISION  DX( 1 ), DY( 1 )
  12. *     ..
  13. *     .. Local Scalars ..
  14.       INTEGER           I, IX, IY, M, MP1
  15.       DOUBLE PRECISION  DTEMP
  16. *     ..
  17. *     .. Intrinsic Functions ..
  18.       INTRINSIC         MOD
  19. *     ..
  20. *     .. Executable Statements ..
  21. *
  22.       IF( N.LE.0 )
  23.      $   RETURN
  24.       IF( INCX.EQ.1 .AND. INCY.EQ.1 )
  25.      $   GO TO 20
  26. *
  27. *        code for unequal increments or equal increments
  28. *          not equal to 1
  29. *
  30.       IX = 1
  31.       IY = 1
  32.       IF( INCX.LT.0 )
  33.      $   IX = ( -N+1 )*INCX + 1
  34.       IF( INCY.LT.0 )
  35.      $   IY = ( -N+1 )*INCY + 1
  36.       DO 10 I = 1, N
  37.          DTEMP = DX( IX )
  38.          DX( IX ) = DY( IY )
  39.          DY( IY ) = DTEMP
  40.          IX = IX + INCX
  41.          IY = IY + INCY
  42.    10 CONTINUE
  43.       RETURN
  44. *
  45. *        code for both increments equal to 1
  46. *
  47. *
  48. *        clean-up loop
  49. *
  50.    20 M = MOD( N, 3 )
  51.       IF( M.EQ.0 )
  52.      $   GO TO 40
  53.       DO 30 I = 1, M
  54.          DTEMP = DX( I )
  55.          DX( I ) = DY( I )
  56.          DY( I ) = DTEMP
  57.    30 CONTINUE
  58.       IF( N.LT.3 )
  59.      $   RETURN
  60.    40 MP1 = M + 1
  61.       DO 50 I = MP1, N, 3
  62.          DTEMP = DX( I )
  63.          DX( I ) = DY( I )
  64.          DY( I ) = DTEMP
  65.          DTEMP = DX( I+1 )
  66.          DX( I+1 ) = DY( I+1 )
  67.          DY( I+1 ) = DTEMP
  68.          DTEMP = DX( I+2 )
  69.          DX( I+2 ) = DY( I+2 )
  70.          DY( I+2 ) = DTEMP
  71.    50 CONTINUE
  72.       RETURN
  73.       END
  74.